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Abstract 

In this paper we analyze the signal emitted by a compact binary system in the 
Jordan-Brans-Dicke theory. We compute the scalar and tensor components of the 
power radiated by the source and study the scalar waveform. Eventually we consider 
the detectability of the scalar component of the radiation by interferometers and 
resonant-mass detectors. 
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1 Introduction 



The detection of gravitational waves (GW) is a field of active research from the point 
of view both of the development of suitable detectors and of the study of possible 
sources and signal analysis. The detectors now operating as GW observatories are 
of the resonant-mass type and have a sensitivity to typical millisecond GW bursts 
of h ~ 6 x 10~ 19 (h is the wave amplitude) or, in spectral units, 10~ 21 (Hz) -1 / 2 
over a bandwidth of a few Hz around 1 kHz |IJ. The first bound is appropriate for 
describing the sensitivity to gravitational collapses while the square of the second 
bound represents the input GW spectrum that would produce a signal equal to the 
noise spectrum actually observed at the output of the detector. With this sensitivity 
it is possible to monitor the strongest potential sources of GW in our galaxy and in 
the local group (distances of ~ 1 Mpc). In order to improve the sensitivity of these 
instruments, more advanced transducers and amplifiers are under development as 
well as new resonant-mass detectors of spherical shape. Furthermore a huge effort 
is under way to build large laser interferometers. It is widely believed that in the 
near future, sensitivities of the order of 10 -23 (Hz)" 1 / 2 over a bandwidth of several 
hundred Hz will be attained allowing the observation of GW sources up to distances 
of the order of 100 Mpc 0. It thus seems that the detection of GWs is highly 
probable at the beginning of the new millennium. In addition to information of 
astrophysical interest, the detection of GW gives an opportunity to test the content 
of the theory of gravity. In fact, it has been shown that a single spherical resonant- 
mass detector , or an array of interferometers |1 , have the capability to probe the 
spin content of the incoming GW. 

One of the most intensively studied GW source is the inspiralling compact binary 
system made of neutron stars or black holes. In the Newtonian regime, the system 
has a clean analytic behaviour and emits a wave-form of increasing amplitude and 
frequency that can sweep up to the kHz range of frequencies. In this paper we 
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study the radiation emitted by this source in the framework of the Jordan-Brans- 
Dicke (JBD) theory. We consider this theory to be of particular interest, since the 
coupling between the scalar field and the metric has the same form of that of string 
theory, which is widely believed to give a consistent quantum extension of classical 
gravity Our main motivation then comes from the attemp to explore a possible 
experimental signature of string theory as already discussed in ||. Furthermore the 
results obtained here generalize to any theory with a JBD type coupling between 
matter and gravitation. 

There has been much work in this domain in the past years. Before going to the 
plan of the paper, we shortly review it. In || binary systems were first proposed as 
possible sources from which extract more stringent bounds on u BD (see (2.1) for its 
definition) than those obtained from solar system data. An analysys of spherically 
symmetric collapse of inhomogeneous dust was carried on in || and later confirmed 
in |J . The case of homogeneous dust was treated in . In || a test particle around 
a Kerr black hole was studied and results very similar to those of our section 5.1 for 
interferometers were found. In || it was pointed out that deviations from general 
relativity can be much different in strong and weak gravity. In these deviations 
were parametrized in a two dimensional space and exclusion plots were drawn out 
of the available data. Finally in |T| spherical collapses were studied in a formalism 
which kept in account strong gravity effects. 

The plan of the paper is the following: in Section two we describe the scalar and 
tensor GW solutions of the JBD theory. In Section three we compute the power 
emitted in tensor and scalar GW by a binary system. In Section four we concentrate 
on the scalar waveform. In Section five, we study the interaction between the scalar 
waveform and two types of earth-based detectors: interferometers and spherical 
resonant-mass detectors, giving limits for the detectability of the signals coming 
from typical binary sources. Eventually, in Section six, we draw some conclusions. 



2 



2 Scalar and Tensor GWs in the JBD Theory 



In the Jordan-Fierz frame, in which the scalar field mixes with the metric but 
decouples from matter, the action reads fT2|| 



Jd 4 X^ QR-^g^d^d^ +~ c J d 4 xL m [^ mjgilu ] , (2.1) 



c 3 



16tt 

where uj bd is a dimensionless constant, whose lower bound is fixed to be uj bd « 600 
by experimental data [ITS], g^ u is the metric tensor, is a scalar field, and ip m 
collectively denotes the matter fields of the theory. 

As a preliminary analysis, we perform a weak field approximation around the 
background given by a Minkowskian metric and a constant expectation value for 
the scalar field 

if = y?o + £ • (2.2) 

The standard parametrization ip = 2(u BD + 2)/G(2u BD + 3), with G the Newton 
constant, reproduces GR in the limit uj bd — > oo, which implies (fo — *■ 1/G. Defining 
the new field 

1 £ 

0^ u = - -rj^h - r)fj, u — (2.3) 

2 y? 

where h is the trace of the fluctuation h^, and choosing the gauge 

d^ u = (2.4) 

one can write the field equations in the following form 

dad^ = -—v (2.5) 
d a d a i = — S (2.6) 
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where 



V, = — (T^ + V) (2.7) 



2 (^o 



(2.8) 



2(2w BD +3)V 2 (f J 16tt 
In the equation ( |2.7| ), T Mi , is the matter stress-energy tensor and t^ u is the gravita- 
tional stress-energy pseudo-tensor, that is a function of quadratic order in the weak 
gravitational fields 9^ u and £. The reason why we have written the field equations 
at the quadratic order in 9^ v and £ is that in this way, as we will see later, the 
expressions for 9^ IU and £ include all the terms of order (v/c) 2 , where v is the typical 
velocity of the source (Newtonian approximation). 

Let us now compute r 00 and S at the order (v/c) 2 . Introducing the Newtonian 
potential U produced by the rest-mass density p 

U(x,t) = [ P f ,t } l d 3 x' (2.9) 
J | x — x I 

the total pressure p and the specific energy density IT (that is the ratio of energy 
density to rest-mass density) we get (for a more detailed derivation, see flip : 



r 00 = -P , (2.10) 

T ( I £ \ 

2(2u BD +3)\ 2 <p Q ) 

l! (l + n-3 P - + 2uJ °° + 1 If) (2.11) 



2(2^ BC +3V P ^ BD + 2 

Far from the source, the equations Q2.5| ) and ( |2.6| ) admit wave-like solutions, 
which are superpositions of terms of the form 

9 flu (x) = A^(x,uj)exp(ik a x a ) + c.c. (2.12) 

£(x) = B(x, uj) exp(ik a x a ) + c.c. (2-13) 



Without affecting the gauge condition (|2.4|) , one can impose h = —2£/(p (so that 
dfiv = h^). Gauging away the superflous components, one can write the amplitude 



A^y in terms of the three degrees of freedom corrisponding to states with helicities 



±2 and |15|. For a wave travelling in the ^-direction, one thus obtains 
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(2.14) 



where b = B /ip Q . 



3 Power emitted in GWs 



The power emitted by a source in GWs depends on the stress-energy pseudo-tensor 
t^ u according to the following expression 



P GW = r 2 J <M = r 2 J <t ok > x k dn 



(3.1) 



where r is the radius of a sphere which contains the source, Q is the solid angle, $ 
is the energy flux and the symbol < ... > implies an average over a region of size 
much larger than the wavelength of the GW. At the quadratic order in the weak 
fields we find 



< t 0z >= -z 



32tt 



4K C + 1) 



< (3)0(3)0 > + < {d h afi ){d h af> ) > 



(3.2) 



Substituting ( ^T2|) , (pl3|) into Q, one gets 
< t 0z >= 



A 2 

, tp c 0J 
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2{2uJBD+3) -\B\ 2 +A^*A a ^- - | A" 2 
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(3.3) 



and using (|2.14|) 

< t 0z >-- 



A 2 

8tt 



en | 2 + | e 12 | 2 +(2cj sd + 3) | b 



(3.4) 



From (|3.4| ) we see that the purely scalar contribution, associated to b, and the 
traceless tensorial contribution, associated to e^ u , are completely decoupled and can 
thus be treated independently. 



3.1 Power emitted in tensor GWs 



Eq. (|2.7|) differs from the corresponding tensor field in GR only by a multiplicative 
factor. Then we can directly write the final result using the well-known expression 
for the power emitted by a system of binary stars in GR 

(Pten)n = -^—{Pqr)* (3.5) 

If we take u BD = 600 in (|3.5|), the multiplicative factor 1/Gip differs from one for 
one part in 10 3 . 

In ( |3.5| ) (Pgr)™, is the power emitted at frequency nu (where ujq is the orbital 

frequency) by a system of binary stars according to GR |TB|, [T7| , averaged over one 

period of the elliptical motion and calculated in the Newtonian approximation 

32 G 4 m 2 m%m . . 
(Pgr)u = -T — ^i^g(^ e) (3.6) 

DC (X 

where mi and m,2 are the masses of the two stars, m is the total mass m = mi + m 2 , 
a is the major semi-axis and e is the eccentrity of the ellipse. The function g(n; e) 
depends on the Bessel functions Jfc(ne) 



i 

g{n-e) = -[J 2 n (ne)(e 2 - 2) 2 /(n 2 e 4 ) + 4£(ne)(l - e 2 ) 3 /e 4 



+ J 2 (ne)/(3n 2 ) + AJ n (ne)J' n {ne){e 2 - 2)(1 - e 2 )/(ne 3 ) 
- 8J n (ne)J' n (ne)(l-e 2 ) 2 /(ne 3 )} 

+ 4J' r 2 (ne)(l-e 2 ) 2 /e 2 + 4J> 2 (ne)(l-e 2 )/(n 2 e 2 ) (3.7) 

In the last phases of the binary system evolution, the orbit becomes more and more 
circular, because the bodies radiate the most at their closest approach ]T6| . In the 
case of null eccentricity e = 0, the function g(n;e) reduces to a Kronecker delta, 
g(n; e = 0) = 5 n 2, and the tensor GW frequency is twice the orbital frequency. 



Summing over all the harmonics n [16|, one obtains 



where 

m - jr^yn (> + T/ + ^) < 3 - 9 ) 

Eqs. (|3.6j ) and ( p.8|) are obtained in the approximation of point-like masses (weak 
self-gravity). For compact binary systems like PSR 1913 + 16, they can be used 



upon replacing the masses, mi,m2, by the Schwarzschild masses of the stars |L8 

3.2 Power emitted in scalar GWs 

We now rewrite the scalar wave solution (|2.13|) in the following way 



t) = £(x, u)e- iujt + c.c. (3. 10) 

In vacuo, the spatial part of the previous solution (|3.10| ) satisfies the Helmholtz 



equation 

(\7 2 + lu 2 )£(x,lu) =0 (3.11) 
The solution of ( 3.11| ) can be written as 

£(x,u) = ^X jm h?\ur)Y jm (6,<p) (3.12) 



where hf (x) are the spherical Hankel functions of the first kind, r is the distance of 
the source from the observer, Yj m (8,(p) are the scalar spherical harmonics and the 
coefficients Xj m give the amplitudes of the various multipoles which are present in 



the scalar radiation field. Solving the inhomogeneous wave equation (2.6), we find 



X jm = WTriuj / j l (ur')Y l l(e,i P )S(x,u) dV (3.13) 



where ji(x) are the spherical Bessel functions and r' is a radial coordinate which 
assumes its values in the volume V occupied by the source. 



Substituting ( p.2|) in (ft.l|) , considering the expressions (|3.10| ) and ( P-12j ), and 
averaging over time, one finally obtains 

(2u BD + 3)c 4 



seal 



E I X jm T (3-14) 



To compute the power radiated in scalar GWs, one has to determine the coefficients 
Xj m , defined in ( |3.13j ). The detailed calculations can be found in the appendix A, 
while here we only give the final results. Introducing the reduced mass of the binary 
system \i = mim 2 /m and the gravitational self-energy for the body a (with a = 1, 2) 

1 r p{x)p{x') « , 3 , 



a ' = -2kW=^ dXdX <3 - 15) 
one can write the Fourier components with frequency uujq in the Newtonian approx- 
imation as (see flQ5|) , ( |OID , flQ3) ) 

(A 00 )n = 5 — nJ n {ne) (3.16) 
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(X 



l±l)n 



{X 2 ±2)n 



3 uj bd + 2 \m 2 mi) 
1 



±J'„(ne 



(l-e^Une) 



2 /IT tLUo^ipo 



3V 5u BD + 2 



\ia nJ n (ne) 



(3.17) 
(3.18) 



T2 



7T IUJ "ip 2 

ua 

30a; BD + 2^ 



n 



-{(e 2 - 2) J n (ne)/(ne 2 ) + 2(1 - e 2 ) J' n (ne)/e 



T 2(l-e 2 ) 1 / 2 [(l-e 2 )J n (ne)/e 2 - 4(ne)/(ne)]} (3.19) 
Substituting these expressions in (|3.14|), leads to the power radiated in scalar GWs 



in the n-th harmonic 

(Pscal)n = Pt° + Ft" + 

where the monopole, dipole and quadrupole terms are respectively 
pj=o . 



(3.20) 



pj'=i 
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m 3 /j, 2 G 4 


n 2 J 2 {ne) 


9K D + 2) 


a 5 c 5 
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m 3 /j, 2 G 4: 


m(n; e) 


9(u BD + 2) 


a 5 c 5 
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m 2 \i 2 G 3 




3K. + 2) 


a 4 c 3 


Vm2 mi 


n 2 


j£(ne) 


1 , 

+ - I" 

e z 


-e 2 )J 2 (ne) 
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m 2 \i 2 G 3 


/fi 2 fii 



(3.21) 



pj=2 



3(lu bd + 2) a 4 c 3 Vm 2 mi/ 

8 m 3 /i 2 G 4 , , 
g{n;e) 



d(n; e) 



(3.22) 
(3.23) 



15(u BD + 2) a 5 c 5 

In figures 1, 2, 3 we plot the monopole m(n; e), dipole d(n; e) and quadrupole g(n; e 
functions against the index n, for different values of the eccentricity e. 
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e=0.7 




Figure 3 
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From the figures one can infer the dominant harmonics in the scalar GWs radia- 
tion. In the case of circular orbit, the dipole function d(n; e) reduces to a Kronecker 
delta d(n; e = 0) = S n i, while the monopole function m(n; e) goes to zero. 

The total power radiated in scalar GWs by a binary system is the sum of three 
terms 

Pscal = P 3= ° + P j=1 + P 3=2 (3.24) 



where 



P j =° = 1 16 X™ 21 "?" 1 , e \, 7 J l + ~) (3-25) 
9{u BD + 2) c 5 a 5 (1 - e 2 ) 7 / 2 V 4 / 



pj=1 2 /Q 2 G 3 m\m\ 1 



e 



■r 



uo BD + 2\m 2 mj c 3 a 4 {I - e 2 f/ 2 \ + 2 ) ^' 2 ^ 
pJ=2 _ 8 G' mjmlm 1 / 73 37 A 

15K D +2)c 5 as (1-e 2 ) 7 / 2 I 24 + 96 7 1 J 

Note that P J=0 , P J=1 , P J=2 all go to zero in the limit u BD — > oo. 

4 Scalar GWs 

We now give the explicit form of the scalar GWs radiated by a binary system. To 
this end, note that the major semi-axis, a, is related to the total energy, E, of the 
system through the following equation 

Gm x m 2 
a = -^2E- 

Let us consider the case of a circular orbit, remembering that in the last phase of 
evolution of a binary system this condition is usually satisfied. Furthermore we will 
also assume mi = m 2 . With these positions only the quadrupole term, ( |3.23| ), of 



the gravitational radiation is different from zero. The total power radiated in GWs, 
averaged over time, is then given by ( |3.8j ), ( p.25| )-( [3~27j ) 

8 G 4 m\m 2 l m 



lh{uJ BD + 2) c 5 d 5 
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■[6(2^ SD + 3) + l] (4.2) 



where d is the relative distance between the two stars. The time variation of d in 
one orbital period is 

d = -9H^p (4.3) 

2E 2 v ; 

Finally, substituting (|4. 1|) , (|4.2|) in (|4.3|) and integrating over time, one obtains 



(2 12^, D + 19G 3 m 1 m 2 m V /4 j4 

rf = 2 (j^^ — ~J (44) 

where we have defined t = t c — t, t c being the time of the collapse between the two 
bodies. 



From ( |3.12| ), ( |3.16| )-( |3TT9| ) and ( [B.3Q one can deduce the form of the scalar field 



(see appendix B for details) which, for equal masses, is 

2/x 



v 2 + - {n-vf + ^(n- d) 
d cr 



(4.5) 



r(2u BD + 3) 

where r is the distance of the source from the observer, and n is the versor of the 
line of sight from the observer to the binary system center of mass. Indicating with 
7 the inclination angle, that is the angle between the orbital plane and the reference 
plane (defined to be a plane perpendicular to the line of sight), and with ip the true 
anomaly, that is the angle between d and the x-axis in the orbital plane x-y, yields 
n- d = d sin 7 sinip. Then from (|4.5| ) one obtains 



m = , 9 2G fTw,/ sin2 7 cos(2V(t)) (4.6) 
which can also be written as 

e(r)=e (r)sin( X (r) + x) (4.7) 

where x is an arbitrary phase and the amplitude £o( r ) is given by 

2G\im . 2 
£0 (r) = j- -——sin 7 

2(2cu flI) +3)r ^12^ + 19^ l2c n y r 1 / 4 r v ; 
In the last expression, we have introduced the definition of the chirp mass M c = 
(m 1 m 2 ) 3 / 5 /m 1 / 5 . 
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5 Detectability of the scalar GWs 



Let us now study the interaction of the scalar GWs with two types of GW detectors. 

As usual, we characterize the sensitivity of the detector by the spectral density of 
strain Sh(f) [Hz] -1 . The optimum performance of a detector is obtained by filtering 
the output with a filter matched to the signal. The energy signal-to-noise ratio SNR 
of the filter output is given by the well-known formula: 

r+oo \H(f)\ 2 

SNR= l -^fHdf (5.1) 

where, in our case, H(f) is the Fourier transform of the scalar gravitational waveform 
h s (t) = GCo(t). 

We must now take into account the astrophysical restrictions on the validity of 
the waveform (|4.7|) which is obtained in the Newtonian approximation for point-like 
masses. In the following, we will take the point of view that this approximation 
breaks down when there are five cycles remaining to collapse [p4 , pB| . 

The five-cycles limit will be used to restrict the range of M c over which our 
analysis will be performed. From (|4.4|) , one can obtain 

u 3 (t) = 2uj = 2\ 



'Grn 
~~dF~ 



(^L\ m ( -» +2 (5.2) 



V64G 5 / 3 / \12lu bd + 19J M c 5/s 
Integrating (|5.2j) yields the amount of phase until coalescence 



Setting ( |5.3| ) equal to the limit period, T 5 cyc i es = 5(2n), solving for r and using ( |5.2 ) 
leads to 

3 /5 

^5 cycles = 2vr(6870 Hz) ( " BD + 2 ) ^ (5.4) 
5 y v ; V12cj bd + 19; M c K ' 
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Taking u BD = 600, the previous limit reads 



^ 5cs/ des = 27r(1547 Hz)^ (5.5) 



5.1 Interferometers 

An interferometric detector measures the relative displacements between the mir- 
rored faces of test masses arranged in the L-shaped configuration of a Michelson 
interferometer. The directivity antenna pattern for a tensorial wave is such that 
the maximum detector output is obtained for a wave inpinging perpendicularly re- 
spect to the plane defined by the interferometer arms. On the contrary a scalar 
BD wave (which is also transverse) inpinging in the same direction will give a null 
effect. In the case of a scalar wave, the maximum effect will be obtained for a wave 
propagating along one interferometer arm. Assuming such a direction and setting 
sin 7 = 1 in ( |4.8| ), we can, for instance, evaluate the SNR for the VIRGO interfer- 
ometer, presently under construction. We use for Sh{f) the VIRGO noise spectrum 
as modeled in [HJ], which is the sum of three main components: thermal noise in 
the pendola, thermal noise in the mirrors, shot noise at high frequency: 



S h (f) = 10 



-47 



( f Y\ ( f Y\ ( f 



100Hz / V 100Hz / V 100H Z/ 

where «i = 2.0, a 2 = 91.8 and a 3 = 1.23. 



Hz" 1 (5.6) 



Integrating ( |5.1|) over the range 10 Hz-500 Hz, one obtains the following SNR: 



SNR = 7.7 10 4 4^ 



2 'M.\ 5/3 / to +2 V /2r 



Mpc/ \MqI V 12w sd + 19 



2(2^ SD +3) 



(5.7) 



For uj bd = 600, we find: 



-2/ M x5/3 



SNR = 3.8 10" 3 — - \^\ (5.8) 
\M VC ) \M & ) 

The inspiralling of two neutron stars of 1.4 solar masses each will then give SNR = 1 
at a source distance r ~ 70 kpc. The inspiralling of two black holes of 10 solar masses 
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each will give SNR = 1 at a distance r ~ 300 kpc. To get this last limit we have 
integrated ( |5.1| ) up to cyc i es = 27r(178 Hz). 



5.2 Spherical detectors 

A GW excites those vibrational modes of a resonant body having the proper sim- 
metry. Most people consider the next generation of resonant-mass detectors will be 
of spherical shape. In the framework of the JBD theory the spheroidal modes with 
/ = 2 and I = are sensitive to the incoming GW |2|, [2(| . Thanks to its multimode 
nature, a single sphere is capable of detecting GW's from all directions and polariza- 
tions. We evaluate the SNR of a resonant-mass detector of spherical shape for its 
quadrupole mode with m = and its monopole mode. In a resonant-mass detector, 
Sh(f) is a resonant curve and can be characterized by its value at resonance Sh(f n ) 



and by its half height width [21]. Sh(f n ) can thus be written as 



SM = %V§T (5 ' 9) 

u n y °cnjn 



Here a n is the cross-section associated with the n-th resonant mode, T is the 
thermodynamic temperature of the detector and Q n is the quality factor of the 
mode. 

The half height width of Sh(f) gives the bandwidth of the resonant mode 



Mn = j^T-^ (5.10) 

Here, T n is the ratio of the wideband noise in the n-th resonance bandwidth to 
the narrowband noise. 

From the resonant-mass detector viewpoint, the chirp signal can be treated as a 
transient GW, depositing energy in a time-scale short with respect to the detector 
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damping time. We can then consider constant the Fourier transform of the waveform 



within the band of the detector and write 21 



SNR 



2nAf n \H(f n ) 

Sh(fn) 



(5.11) 



The cross-sections associated to the vibrational modes with I = and I = 2, 
m = are respectively f22| 



C(rcO) 
C(n2) 



H, 



GMv< 



V( Wm +2) 
F n GMv? 



(5.12) 
(5.13) 



6 d»(w m +2) 

All parameters entering the previous equation refer to the detector M is its mass 



v s the sound velocity and the constants H n and F n are given in |[22|| . The signal-to 
noise ratio can be calculated analytically by approximating the waveform with a 



truncated Taylor expansion around t — 0, where a> g (£ = 0) = c<j n ; [24 



/i s (t) w G£ (t = 0) sin 



1 (du\ 



(5.14) 



Using quantum limited readout systems, one finally obtains 

5 • 2 1 / 3 H„G 5/3 



(SNR 



n)l=0 



32K D +2)(12^ flD + 19)fcc3 



M r 5/3 Mv* 



r 2 ^ n o 4/3 



■ 4 

sin 7 



(5.15) 



(SNR, 



5 • 2 1 / 3 F r ,G 5 / 3 



n 1=2 



192(uo BD + 2)(12uj BD + 19)hc 3 

M c 5 / 3 Mv s 2 . 4 
■ sin 7 



(5.16) 



which are respectively the signal-to-noise ratio for the modes with I = and I = 2, 
m = of a spherical detector. 

It has been proposed to realize spherical detectors with 3 meters diameter, made 



of copper alloys, with mass of the order of 100 tons [2c]. This proposed detector 
has resonant frequencies of u>u = 2n ■ 807 rad/s and u>\q = 2n ■ 1655 rad/s. In the 
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case of optimally oriented orbits (inclination angle 7 = tt/2) and uo BD = 600, the 
inspiralling of two compact objects of 1.4 solar masses each will then be detected 
with SNR = 1 up to a source distance r(cj 10 ) ~ 30 kpc and r(co> 12 ) — 30 kpc. 



6 Summary and conclusions 

In this paper we have studied the waveforms emitted by a system of binary stars in 
the framework of the JBD theory and computed the power emitted in GW's for the 
tensor and scalar components. Eventually we derived limits for the detectability of 
such signals by interferometers and resonant mass detectors. In the former case we 
left aside the question of the detectability of the scalar component of the GW || 
and we have concentrated on waves impinging from the most favourable direction. 
We would now like to comment The detectability ranges obtained in Sections 5.1,5.2 
for the scalar component of the GW's emitted by a binary system, vary from few 
tens to few hundreds of kpc's for masses ranging from those of typical neutron stars 
(1.4M ) to those of typical black holes (1OM ). We remind the reader that for the 
purely tensorial component (in this case the results obtained in the framework of 
the JBD theory are practically the same of those of GR) the detectability range 
(for 1.4M Q ) is r ~ 120 Mpc for spherical detectors and r ~ 300 Mpc for 
interferometers ||19|| . The expected rate of coalescence events is of the order of 1 



per year up to 100 Mpc p7 |. We can thus conclude that binary systems look less 



promising than gravitational collapses |22] as sources of detectable scalar GW from 



the next generation of earth-based detectors. 
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A Appendix 



In order to calculate the coefficients Xj m defined in ( [3.13| ), let us first express Xj m 
as a sum of Fourier components 

Xjm(t) = E (X im )ne inf3 (A.l) 

n=— oo 

where (3 is the mean anomaly which, in terms of the orbital frequency ujq and of the 
time of periastron passage T (or equivalently in terms of the eccentric anomaly a 
and the eccentricity), results 

(3 = u)o(t — T ) = a — e sin a (A. 2) 

In the so-called quadrupole approximation y = uod/c <C 1, the spherical Bessel 
functions ji(y), which appear in (|3.13| ), can be written as 



Jliy) = 2 2*fc!(2Z + 2fc + l)!! (A ' 3) 

Making use of the Newtonian approximation (including only two terms of the series 
( |A.3[ )), from ( J2.ll ) and ( |3.13j ) one can obtain X o as 



u} BD +2Jv\ 6 J \ p u BD + 2 J 



4) 



To simplify ( |A.4 ) we use the post-Newtonian expressions of the conserved quantities 




p° = fp(i + v 2 + 5uJbd + 4 N u + n )dx 

J v 2 k c + 2) ) 



(A.5) 



.2 , 5UJ BD + 4 T T i TT i P 



p l = p[i + v z + - BD - u + n + r -\ v l dx 
1 V 2 K D + 2) pI 



- \fp(l + t + 3± ^U)W*dx (A.6) 



2 u BD + 2 



where 





X 


— x! \](x — x')i 




x — x ' 


3 



W* = / p' [U ' 1 ;i dx' (A.7) 
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To the required accuracy and modulo constants, one then obtains 



y/4mu<p fj, 2 (udf m\ 
u BD +2 \ 6 rj 

In terms of the eccentric anomaly a, the above expression reads 



\ B D + 2 ) 

where 

l + ecos« n 2 2 1 

G i=0 =- + — (1 - ecosa) 2 + (A.10 

1 — e cos a o 1 — e cos a 

Let us express Gj=o as a sum of Fourier components 

oo 

G j=0 = J2 A n cos(np) (A. 11) 

n=0 



where (3 is the mean anomaly defined in (|A.2|) . Using the integral expression of the 
Bessel functions 

1 r n 

Jn{%) = — / cos(na — x sin a)da (A. 12) 

7T Jo 



and the recursion formula 



2n 

J n +i(x) + J n -i(x) = — J n (x) (A. 13) 

X 

results 

2 /•"" 16 
A n = - G j=0 cos(n(3)d(3 = — J n (ne) (A.14) 

7T JO O 

Then one finally has 

^v 7 ^ zw ^o m// 

(Aoo)n = ^ nJ n {ne) (A. 15) 

3 u; SD + 2 a 

Choosing the orbital plane as the x-y plane, yields X 1Q = and 

^-^/^W i+n -^^) (A - i6) 

Defining Q = 0,2 /m^ — Oi/mi, one then obtains to the required order 

X 1±1 = ^f^ g d( C o S ^T^sin^ (A.17) 
VQrr{u BD + 2) 
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where ip is the true anomaly. In terms of the eccentric anomaly a, (|A.17 ) results 

Xl±1 = (A.18) 

V 6 BD + 1 

where 

G j=1 = ±(cosa-e) - z(l -e 2 ) 1/2 sina (A. 19) 

The binary system center of mass calculated with respect to the gravitational self- 
energies Q a doesn't coincide with the center of mass with respect to the inertial 
masses m a of the two bodies, if the masses are different (Nordtvedt effect): the 
resulting dipole moment is, as we have seen, a source of scalar radiation. If we 
express Gj=\ as a sum of Fourier components 

oo 

G j=1 = B + i B n cos(np) + C n sm(nP)] (A.20) 

n=l 

we obtain 

B n = ±-J' n {ne) (A.21) 
n 

(where the prime indicates derivative with respect to the argument ne) and 

C n = - £ G J=1 mi(n(3)d(3 = -^(1 - e 2 ) 1/2 J n {ne) (A.22) 
Then the n-th component of the coefficient X\±\ is 

[X 1± i) n = -\ — rr;yp a 



3 u BD + 2 
±J' n {ne)--{l-e 2 ) 1/2 J n {m i 



e 

Finally, in the case j = 2 one obtains X 2 ±i = and 



(A.23) 



X.o = TTp~~7y{ i r' 2 Y^9A) P dV (A.24) 
lh{u BD + 2) Jv 

X 2±2 = 4 ^Vo / r ' 2 Y* ±2 {eA)pdV (A.25) 
15 co> + 2 Jv 



lh(co BD + 2) Jv 
and in terms of a, using the Newtonian approximation 



x w (t) = \ -ir ( -^(i-cos") (A- 26 ) 

X 2±2 = T\ ^ Wa W 2 (A.27) 

* V30(cj sd +2) J " 2 v ; 
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n=l 



where 

oo 

G j=2 = D + J2 i D n cos(n(3) + E n sm(n{3)} 
Calculating D n and E, 

D n = 

n \_ne z v 

1 



4 



' (2-e 2 )J n (ne) + -(e 2 -l)J' n {ne) 

e 



n 



ne 

"1-e 2 



3 -Jn{ne) 



tie 



-J' n (ne) 



the n-th components result 



2 fW iu> Q 3 (p 2 



3V5o; SD + 2 



/xa nJ n (ne) 



( v s „ /7T HTO 2 

(^2±2jn = =F2\/ — —//a 



30^ flD + 2' 



n 



-{(e 2 - 2)J n (ne)/(ne 2 ) + 2(1 - e 2 ) J>e)/e 



B Appendix 



(A.28) 

(A.29) 
(A.30) 

(A.31) 



T 2(l-e 2 ) 1 / 2 [(l-e 2 )J n (ne)/e 2 - j;(ne)/(ne)]} (A.32) 



We want to determine explicitly the form of the scalar GWs radiated by a binary 
system. From ( [Q2]) and QPP , (|Q7| ), ( TA~2^ ), and taking into account that 

in the limit r — > oo the spherical Hankel functions become 



3 i[ajr-(/+l)7r/2] 



(B.1) 



one can easily obtain 



2/i 



(2u; flD +3)Gr 
{v 2 + m/d + (cud) 2 /6 + 2ito(^l2/m2 — £l\/mi)n ■ d 

-u 2 d 2 (3n 2 - 1)/12 + cu 2 [(n • cf) 2 - (n,d, - 7vk) 2 ]/4} w (B.2) 

where the subscript uo indicates that all the quantities in the right member of the 
above expression ( |B.2| ) are to be considered as Fourier components with frequency 
oj (for example, v — ■> u(tt')). 
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The time dependent amplitude is [13 



BD 

2 m 1 rf 2 , , . ks /„ 3n z 2 - T 

u + -r T^(dkdr) 1 ! 

d 6 dt 2 v ; I 2 



(2o, SD +3)Gr 



d 2 



2£n-t;--^(n-d) 2 



[n x dy fiyd x ^ 



2\i 



(2u BD +3)Gr 



f H — - — 2£n • t> — (n • t> ) + — n ■ a 



d 



d? 



(B-3) 



where we have set v = v(uj)e tu} ^ ^ + c.c, ecc. 
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